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I. INTRODUCTION 

The customary procedure to deal with charged particle scattering is to partial wave the 
amplitude and identify the Coulomb phase shift from which scattering information can be 
obtained, by adding to the Coulomb amplitude, the contribution from whatever other short- 
range potential. In many applications, the asymptotic form of large angular momentum 
is employed for the Coulomb phase shift. In this paper we revisit the derivation of this 
asymptotic and derive next to leading order correction to the usual WKB from. 
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First, we recall the form of the Coulomb phase shift pQ 

2lai _ Tjl + l + ty) 

T(l + l-ir,y [) 

where T is Euler's gamma function, /, a non-negative integer, is the angular momentum and 
the real parameter 77 is the so-called Sommerfeld parameter, which is inversely proportional 
to the square root of the scattering energy. 

In this paper we will present simple formal proofs of some asymptotic approximations for 
the gamma function, like Stirling series and Gudermann series, and we will discuss several 
methods for approximately computing the phase shift, Eq. (fij), from these asymptotic 
approximations in various regimes. We will compare these results with the corresponding 
ones obtained from other methods for computing phase shifts, like the WKB and the eikonal 
approximations. We will also present proofs of some exact relations for the phase shifts 07 
which are more or less known in the literature. For instance, using Eq. ([!]) we will easily 
show that 

1 

V 

tan _ 1 

m=l 



I 

^ = ^o + ^tan- 1 (^), (2) 



and that 



<?o = -IV ~Y] tan 1 ( j 
' V m + 1 / 



-1 V \ V 

Ifclil 

m=0 

valid for all rj G R, from which we obtain for \rj\ < 1 the power series representation 



(3) 



'.^-E™^. (4) 



2k + 1 
k=i 



where ( is Riemann's zeta function and 7 is Euler's constant. 
Using our asymptotic expansions we show that 



4 1} = \ tan" 1 ( V ) + V (in (^1+^) " l) " n (5) 



is an excellent approximation for the exact expression ([3]). We also show numerically that 
this approximation is very good even for small values of rj. More generally, we show that 

approximates Q very well for + l) 2 + i] 2 "large" . 
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II. ASYMPTOTIC PROPERTIES OF THE GAMMA FUNCTION 

We start with the well-known integral representation for Euler's gamma function, valid 
for Re(z) > -1 0, 

poo 

T{l + z) = / dtt z e- 1 . (7) 



o 



For ([T]) we will take z — I + %r\. Writing t z = e zlnt , we can use the saddle point method to 
evaluate the integral 

POO 

T(l + z)= dt e ~ t+zln{t) . (8) 



o 



Expanding the exponent of the integrand around its extremum, at t = z and keeping terms 
up to second order, 

-t + z\n(t) ~ -z + z\n(z) - ^-(t-z) 2 , (9) 

Eq. (pi) becomes 

POO 

T{l + z)~e~ z z z dte- {t ~ z)2/2z . (10) 
Jo 

For Re(z) = I 3> 1 and lm(z) — r\ > 0, the Gaussian integral can be easily evaluated as 

POO 

dte - [t - z? /2 z = / dte -(t-^/2 z = ^ (u) 



J -oo 

and using this result in the previous equation, we obtain 

r(l + z) ~ V2^e~ z z z+1/2 . (12) 

This is the well-known Stirling's approximation for the gamma function, whose validity 
can be established for the whole complex plane, except the non-positive real axis, i.e., in 
C~ = C \ (— oo, 0], and for "large" values of \z\. See, e.g., Ref. [2]. 



On the real line, Stirling's approximation (12) has a long history having been first pre- 



sented by de Moivre around 1730, who found an approximation for n!, valid for "large" n, 
in the form n\ ~ K n n+ ^e~ n , for some unspecified constant K. In the same year, Stirling 
proved that K = \f2n using Wallis product formula for it. This approximation became 
known as Stirling's approximation for n\ and quite soon generalizations for Euler's gamma 
function on the positive real line became available. Stirling's approximation is very useful 
in Statistics and Probability Theory because if offers a very good approximation for "large" 
factorials. For small values of n, however, correcting factors are in order. Stirling found 
corrections for the approximation that bears his name for n\ in terms of an asymptotic 
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(but not convergent!) series in 1/n that became known as Stirling's series. Further anal- 
ysis of those corrections for Euler's gamma function T(z), valid on the complex half-plane 
Ke(z) > 0, have been performed by Binet in 1839 (Ref. [32]). Another important contri- 
bution was made by Gudermann in 1845 (Ref. [13]), who found another expression for the 
correcting factors in terms of a convergent expansion of another kind. The most important 
contribution to the study of corrections to Stirling's approximation on the complex plane 
was the work of Stieltjes, dated of 1889 (Ref. [H]), who generalized Stirling's approximation 
and Gudermann's corrections to the whole complex plane, excluding the negative real axis 
(i.e., to C~). For a more detailed account of the developments on the complex plane, see 
[2] and [11] . For some recent contributions to the correcting factors to factorials, see [15] 
and [16J. This last reference contains a list of historical results on corrections on Stirling's 
approximation. 

Let us briefly describe the ideas behind Gudermann's corrections and present Stirling's 
series. Since Euler's gamma function satisfies T(z + 1) = zT(z), we get from (12) the 
approximation 

T(z) ~ V2^e- Z z z ~ 1/2 . (13) 



However, ( 13 ) contrasts with the expression obtained from ( 12 ) itself by replacing z by z — 1 : 

T(z) ~ v^e^- 1 ' (z - l) z - 1/2 . (14) 



Since both expressions (13) and (14) are only valid for \z\ very large, there is no practical 
difference between them. Nevertheless, one can better deal with this situation by seeking an 
exact representation for T in the whole region C~ (and not only for "large" \z\) in the form 

e-' z z ~ 1/2 e^ z) , 



r z) 



(15) 



and fixing the correction factor by imposing the relation T(z + 1) = zT(z). A simple 
computation reveals that this condition implies that \x has to satisfy the functional equation 



H(z) -fi(z + l) 



z + 



In 1 + 



1. 



(16) 



Moreover, the validity of (13) for "large" \z\ leads to the condition limui^oo /i(z) = 0. This 



allows to a solution for (16). Indeed, it follows immediately from (16) that for any positive 
integer n one has 



n-l 



H(z) - n(z + n) = Y^ 



m=0 



z + m + 



In 1 + 



z + m 



(17) 
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Hence, the condition lim^i^oo /j,(z) = implies, in particular, that lim r 



we get 



m=0 



: -h ni + - | In ( 1 



1 



z + m 



n{z + n) = and 



(18) 



This series is known as Gudermann series. According to [2], it was first obtained by that 
author in 1845 [13] for real and positive z and the generalization for z G C~ was obtained 
by Stieltjes in 1889 [H]. 



The series ( 18 ) converges for all z G C and fi can be bounded by 

1 



Hz)\ < 



1 



1 



(19) 



12cos 2 (yV2) |*| ' 

with z = \z\e llf , z G C~ (see [2] or [H]). On the real line, very precise upper and lower 



approximants of the T-function can be obtained from (15) and (18), valid also for small 
values of the argument (see [15] and [T6]). 



Notice that, by (19), for \z\ ~ 1, the contribution of the factor e^ z ' to (15) is lower than 



ci. 9% for if ~ and lower than ci. 18% ip ~ tt/2. Hence, for the sake of precision, it is 



relevant in that range to consider corrections to Stirling's approximation (12) 



With (15), and since T(z + 1) = zT(z), we can also write 



r(z + i) = v^ e -*z* +1/ V w . 



(20) 



Hence, acts as a correcting factor for both (12) and (13). Beyond the Gudermann series 



( 18 ), the function /i can be represented in many other forms. One of the most useful of them 



is the so-called Stirling series: 



H{z) 



E 



B 



2n 



1 



' (2n- l)2nz 



2n-l ' 



(21) 



where Bk is the fc-th Bernoulli number. One has to say that the series on the r.h.s. of (21) 



is an asymptotic series in C , but is not convergent! See again [2]. Therefore, it is not to be 
seen as a Laurent expansion for /i around z — 0. In fact, z = is not an isolated singularity 



of /i, but a branch point, as we see from (18). Since (21) is an asymptotic series we can get 



good approximations for /x in C by truncating it at some finite value of the index n and 



taking l/\z\ small enough. The first terms of (21) are 

1 1 



fi(z) 



1 1 

12 i 



360 2 3 



1 1 
h 1260^ "•■ 



(22) 
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For the correcting factor e^ z \ this gives 



1 1 



1 1 



139 1 



12 z 288 z 2 51840 z 3 



+ 



Therefore, we may write 

T(z) = V2^e- Z z z - 1/2 
r(z + l) = V2^e- z z z+1/2 



1 + 
1 + 



1 1 


+ 


1 


1 


139 1 


12^ 


285 


z 2 


51840 z 3 


1 1 


+ 


1 


1 


139 1 


12z 


285 


z 2 


51840 z 3 



+ 



(23) 

(24a) 
(24b) 



These approximations are also known as Stirling's series for the gamma function. They 
are asymptotic (but not convergent!) expansions for T, valid for z G C~ and l/\z\ "small". 



Although (24) are asymptotic but not convergent approximations for T, they can in some 



practical situations be more useful than the representation (15) or (20) with the convergent 



Gudermann expansion (18). 



III. STIRLING'S SERIES. A FORMAL DERIVATION 



The Stirling's series (24) leads to more accurate estimates of the phase-shifts ([T| than 



Stirling's approximation ( 12 ), even when the value of / is not "large" . There are many proofs 



of Eqs. (21) or (24) in the real or in the complex domain (see, e.g., Ref. [2 [3]), but they are 



all rather involved. We will present now a simple derivation of Eq. (24b) (see also jl]). 

Consider the curve C in the complex t-plane parametrized by s G (— oo, oo) and satisfying 
the parabolic functional mapping [H El E] , 

t(s)-zln (t(s)) - A{z) = S -, 

with A(z) = z — z In z and with t(0) = z. Defining u(t) = t — z ln(t), we can write 

u(t(s))=A(z) + S -. 



One can choose t(s) satisfying t(s) — > as t(s) ~ e^ai for s — > — oo and t(s) ~ y + 
ilm(z) In for s — > +oo. By a continuous deformation of the t-integration curve in (7) 
from the positive real axis to C and by carefully extending the integration to infinity, one 
can write, by Cauchy's theorem, 

2/2 dt{s ^ 



T(l + z) 



dt t z e 



-z+z In 2 



C 



ds 



ds. 
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The derivative dt/ds can be written as a power series in s 2 using the mapping equation, 

dt 



ds 



= E 

k=0 



a 2 k s 



2k 



(25) 



It is a simple matter to evaluate the coefficients a^, by repeated differentiation of the mapping 
equation. We calculate below the first two terms, a , and a 2 (higher order terms can be 
derived similarly). Calling u(t(s)) = v(s), we write 



dv(s) 

ds 
d 2 v(s) 

ds 2 
d 3 v(s) 

ds 3 
d A v{s) 

ds 4 



-(}- 



z\ dt 
tJ ds 



s, 



z\ d 2 t z ( dt\ 



+ 



-r =1, 



t J ds 2 t 2 \ds J 

d 3 t 2z d 2 t dt 2z fdt\ 

t) d¥ + ¥ d¥ d~s ~ ¥ \~ds~) 

dH Az d 3 t dt 3z ( d 2 t 



0, 



ds 4 + t 2 ds 3 ds + 



t 2 \ds 2 



+ 



0. 



_12z d 2 t Ut^ 2 
'l^Jd 2 [ds. 

We now evaluate these equations at the extremum point, defined by the condition 

~dw(ty 



6z (dtX 
¥ Vds) 



dt 



0, 



(26) 



to 



which is t — z and s = 0. The above four equations lead to the results 



dt 
ds 



t=t 



dH 
ds 3 



1 



t=t 



and 



dT 
ds 5 



1 



t=to 



36 z 3 / 2 ' 



(27) 



Thus the coefficients a 2 k above are, 



a 



a 2 



1 

2! 



(Z4 



4! 



"eft" 
ds 

dh 

ds 2 

dh 



= \/z 



t=t 



t=t 



ds 5 



1 



t=to 



1 



4! 36 z 3 ' 2 ' 



Using the Gaussian integral formula, 

dse~ s2/2 s n ds = Vto(n- 1)!!, 



/ 



(28) 

(29) 
(30) 

(31) 
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for n even, the full integral in the T function can be written down in the form of the Stirling 
series above, namely, 

oo 

I\z + 1) = V2^e- Z z z ^(2fc-l)!!a 2fc . (32) 



A:=0 

Thus, evaluating the coefficients a>2k and inserting above, we get 



r(z + l) = V2^ e~ z z z+1/2 



1 1 

+ Y2z + 288 z 2 H 



(33) 



The use of the quadratic mapping analysis of integrals and the generation of appropriate 
asymptotic series in general was demonstrated in, e.g., [7] for the case of the Gamow integral 
employed in nuclear astrophysics. This integral has the general form, 

1(a) = J dxS(x) exp (— x 1/2)' (^4) 

where the function S(x) is usually a slowly varying function of x and can be written as a 
sum S(x) = YlkLo Sk% k , which would then result in a series representation of the Gamow 
integral similar to what was done for the T function above. 

We now use the above results to evaluate the Coulomb phase-shifts (TlJ) . 

IV. AN EXACT EXPRESSION FOR THE PHASE SHIFT AND FIRST 

APPROXIMATIONS 

From ([T]) and (15), and using the identity 

hf^U^P), (35) 



x — xy J \x, 
valid for x, y G M with x > 0, we get after some elementary computations the following 
exact expression for the phase shift: 

a^af + M^ (36) 



where we define 



erf + tan- 1 (jl^ + rj (in (V(l + ±) 2 + V 2 ) - l) , (37) 



and Mi^ = ^ (//(l +1 + irf) — /i(l + I — if])). Below, we will use the series expansion (18) to 



find closed expressions for M/ )7? , for ai and, in particular, for cr . Before we proceed let us 



make some comments about some useful approximation we can obtain from (36)-(37). 
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Using (19) with \z\ = \/ (I + l) 2 + f] 2 and (p = tan 1 (j+i), one finds the bound 

1 



M lt „\ < 



e(/ + i + v/(/ + 1)2 + ^ 



Hence, for I ^> 1 or \r]\ 3> 1 the contribution of M; ir) to (36) can be neglected and we can 
restrict clS db first approximation to a; . 



For I = 0, for instance, one gets from (37) the approximation 



ct ^ o"q 0) = 2 tan 1 fa) + ^ ( ln (a/ 1 + »?) ~ 1 



(38) 



with an error bounded by §fl + a/1 + f? 2 ^ . For |t/| ^> 1, Eq. (38) gives the approximation 

ob^ j +77(111(77) -1). (39) 



For M- <S 1 and I > 1 we get from (|37|) the approximation 



01 ~ 77ln(Z + 1). 



(40) 



V. CLOSED EXPRESSION FOR THE PHASE SHIFT 



Now we will try to find closed expressions for 07 and a . According to Q, using the fact 

that T(z + 1) = zT(z) for all z G C~, we have 

e 2 iCTi = T(l + l + i V ) = (l + tr] )...(l + tr ] )T(l + tr ] ) = (I + irj) ■ ■ ■ (1 + zr?) ^ 
r(l + 1 — 17]) (l — irj) • • • (1 — if]) r(l — W7) (/ — zry) • • • (1 — irj) 



Therefore a t = a + ^ Z)L=i m ) and usin S d 35 h one nas 



m=l 



tan" 



1 ^ 



(41) 



Eq. (41) is a remarkable expression, since it shows that o\ differs from cr by a finite sum. 



Let us now find a more explicit expression for a . According to (36)-(37) 



^0 = \ tan" 1 (77) + r] (ln (Vl + 77 2 ) - l) + M 0ir? , 



(42) 



Now we analyze M 0jr? = ^(/-t(l + 2*7) — Ml — ^v)) niore closely. According to (18) (with 
the change of summation variable m — » m — 1), we have, 



Mr 



0,7? 



1 OO 

m=l 



m + irj H — J In (l H 

2 / \ m + irj 



■ 1 , , 1 

m — 27/ H — In H 

2 / \ m — irj j 



(43) 
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After simple rearrangements and using we can write ( 43 ) as 



1 I" ( l \ f, (m + l + ir)\ , 

M „ = - > m + - In - In 

M 2z^ V 2A \m + l-i V ) 



m + irj 
m — ir\ 



+ irj | In ((m + l) 2 + r/ 2 ) — In (m 2 + r/ 2 ^ 

(44) 



Using (35) and denning 

1 / m + it] 

Arr, = — In 



2i \m — irj 



tan 



and 5 m = \ In (m 2 + V 2 ) 



we can write (44) as Mo, v = limTv^oo Mq „, where 

N 



<, = E 



m=l 



+ ^ ) (Aw+i — An) + B m+ i — B m 



Now, 



and 

N 



N 



^3 (5 m +i — B m ) — Bn + i — B\ 



(45) 



(46) 



m=l 



- ( 1\ - 

2^ ( m + - - An) = J] 



m=l 



m=l 



(m + 1) + - An+i - to + - M 



E- 4 

m=l 



m+1 



N + 1 + -) A N+1 



( l\ N 
( 1+ 2 ) A-^An+x- 

^ ' m=l 



(47) 



Hence, collecting the results, we have 



Mr 



0,7? 



lim 

AM-oo 

3 
~2 



iV + - tan" 
2 1 



x 7 m=l x 

tan" 1 (77) - | In (l + r/ 2 ) 

3 r N ( \ 

7] tan" 1 (r)) - - In M + rf \ + lim 77 m(iV + 1) - ^ tan 1 ( — " — ) 

m=l x y 



.(48) 



Now, adding and subtracting rj Y2m=o ~7+T ^° ^he t erms i n brackets, whose limit is being 



taken in (48), we get 

N 

^(iV+lj-^tan" 1 

m=l 



V 



m + 1 



N 1 

ln(jV + l)- V 

v ' ^ m + 1 



m=0 



With this, (48) becomes 



M 0tV = V (l - 7) - - tan" 1 (77) - | In (l + r? 2 ) - ^ 

m=0 



tan 



E 

m=0 



1 / »7 



tan 



-1 



V 



7] 



m + 1 I m + 1 



+tan 



V 



m + 1 J m + 1 



(49) 



11 



a /n 
0.8- 




FIG. 1: The graph of o"o/7r according to the exact formula (50) for < r] < 4. Notice that gq 
vanishes at rj ~ 1.810. 



where 7 = uhia^c 



E 



JV 

m=0 m+1 



1 



-ln(iV + l) 



is Euler's constant 7 ~ 0.5772156649 . 



Inserting (|49|) into (|42|) we finally get, after some trivial cancelations, 

V 



era 



-7V 



E 

m=0 



tan 



m + 1 



?7 

m + 1 



(50) 



By recalling the Taylor expansion of tan 1 about 0, 

(~l) fc 



tan 



2fe+l 



(51) 



we can write (50) for |?7| < 1 as 

r 1)k 



(TO 



-7^ 



EE 



2 A; 



m 



2fc+l 



00 

E 

fc=i 



(_l)fc C ( 2fe + l) 

2k + 1 ' ' 



m=0 fc=l 

where C is Riemann's zeta function. 

In Figure [T] we plot <jq for < rj < 4. As we discuss below, for larger values <jq follows 



very closely the behavior dictated by the asymptotic approximations (38) or (39). In Figure 



[T]we see that <Jq vanishes at 77 = and t] ~ 1.810. This latter value of rj is coincidentally 
quite close to the value of a/2, which was obtained in [T7] in connection with the scattering 
of identical charged particles (Fermions or Bosons). The Mott cross section, at this critical 
value of rj = y/2, was predicted to be isotropic over a broad range of angles around 90°. 
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VI. THE PHASE SHIFT AND STIRLING'S SERIES 



The above asymptotic formulae (37), (38), (39) and (40) will be discussed further below 



within the low-energy, large-??, large-/ case of the WKB approximation, and the high-energy, 
small-??, large-/ case of the eikonal approximation. In this section we use Stirling series for 



the gamma function (24) to further improve those approximations. 



A first order correction to Eq. (|37j) can be easily evaluated using Stirling's series (22)-(24) 
For this purpose, we rewrite Eq. Q in the form 



e 2ia ' ' F(z) 



with z = 1 + / + if] and of 1 given in (|37j), where 



F(z) 



exp 



1 (\ 



1 



1 

360 



1 



y*\3 



(52) 



(53) 



12 \z z* j 3bu \{zy [z 
where the * symbol refers to complex conjugation. The first order approximation for F(z 



is 



F^(z) = 



exp 



1 fl 



12 V z 



--: e 2iAa ' 



Evaluating the above expression, we find 



Ad 



(o) 



12((/ + l) 2 + r/ 2 )' 

Therefore, the first order approximation to the Coulomb phase-shifts is 



(54) 



(55) 



cr 



(i) 



a 



(0) 



/ + - I tan 



-i 



12((/ + l) 2 + r/ 2 ) 



(56) 



Figure [2] shows Coulomb phase-shifts versus /. 



From (56) we get, in particular, the first order correction to (38) due to Stirling's series: 



4 0, + a4°' 



^ tan" 1 (77) + 77 ( In (yjl + rj^ - 1^ 



V 



12 (1 + r/ 2 )' 



(57) 



It is interesting to compare the approximation (|57|) to our exact expression for do given 

ter 



in (3). Figure [3] we plot the relative error a ° ^ — for values of the Sommerfeld parametf 
between and 5. It shows that is an excellent approximation for <jq, with relative errors 
below 1%, even for "small" values of 77, except, perhaps, near 77 ~ 1.81, where a vanishes. 



13 




FIG. 2: Coulomb phase-shifts as functions of I. Results normalized with respect to tt are given for 
two values of the Sommerfeld parameter. 



0.015 - 



0.01 



0.005 - 




FIG. 3: The relative error — for < r] < 5. The sharp peak around r] ~ 1.81 is due the the 
vanishing of o"o at that point. 

VII. DISCUSSION OF THE RESULTS FOR STIRLING'S SERIES 

The lowest orders approximations (in 1/z) of the previous section are supposed to work 
for \z\ ^> 1, which means / 3> 1 and/or 77 ^> 1. They are, in fact, very accurate, even when 
these conditions are not well satisfied. A first illustration of this fact is presented in Figure 
[ij where we compare the exact phase-shifts (solid line) with the lowest order approximation 
Of (open circles) for a small 77, as functions of I. They are very close. The only exception 
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FIG. 4: The approximate Coulomb phase-shifts <r, and oy for r/ = 0.1, as functions of I. For 
details see the text. 

is the case of Z = 0, where they are quite different. 

For larger values of rj, the agreement is much better. This can be seen in Table [TJ, where 
we compare the two approximations of the previous section with the exact Coulomb phase- 
shifts. Note that the is very accurate even for I = and 77 = 0.1. 

This point can be seen more clearly in Figure [5j where the exact s-wave phase-shifts are 
compared with the approximations a^ and 0y= o . The comparison indicates that the usual 
large-?? approximation, cr^ , is very poor. The situation is completely different with the 
improved approximation, cr^ , which includes the first order correction of Stirling's series. 
In this case, one gets accurate results for any value of the Sommerfeld parameter, rj. 



VIII. COMPARISON WITH OTHER APPROXIMATIONS 

We turn now to well known approximations used to calculate the phase shifts, motivated 
by the physical conditions. The first such approximation is the WKB one invoked to consider 
scattering under semi-classical conditions of short local wave lengths. In this approximation, 
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(0) / (±) I 



(i) 



exact 



/7T 



0.1 -0.01581 -0.01844 -0.01825 

1 0.01413 0.01346 0.01348 

2 0.02967 0.02930 0.02938 



1.0 -0.08299 -0.09625 -0.09602 

1 0.1592 0.1539 0.1540 

2 0.3042 0.3015 0.3016 

TABLE I: Coulomb phase-shifts as functions of I, for two values of the Sommerfeld parameter. 



one finds for the phase shift the following expression pQ, 



oi = lim 



dr'ki(r') — 



dr'ki (r') 



(0) 



(58) 



where fcj(r) and kf\r) are the local wave numbers in the presence and absence of the 
potential, respectively, 



Hki(r) = pi(r) 



'2/i (e-V(t) 



mj. + 1/2)- 
2/ir 2 



hk l (r) = p^'ir) 



(0), 



2/i [E- 



h\l + 1/2) 2 

2^ 

(o) 



(59) 



(60) 



where /i, here, is the reduced mass. The radii, r\ and r\ , are the classical turning points 



defined by pi(ri) = 0, and p\ 



i 



0, respectively. The phase shift, now a function of 



the energy (or the asymptotic wave number k), and the semi-classical angular momentum, 
A = Z+l/2, can be evaluated once the potential is given. In the case of point-charge Coulomb 
scattering, the result of such a calculation results in the following expression for the WKB 
phase shift function 5(\, k) = a(X, k), 



tr(A, k) = -r] In [r] 2 + A 2 ] + A sin -1 



-rj. 



(61) 
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FIG. 5: Coulomb s-wave phase-shifts as a function of the Sommerfeld parameter. Exact values are 
compared with the approximations discussed in the text. 

where the Sommerfeld parameter rj is related to the asymptotic wave number, k, by, rj = ka, 
where a is half the distance of closest approach for head-on, zero impact parameter, collision. 



The above expression should be compared to our Eq. (37). Clearly the major difference 



resides in the introduction of the semi-classical angular momentum variable A in the former 



equation Eq. (61 ). This change is in fact necessary when the WKB form of the wave function 
is invoked |Sj, to guarantee the presence of a classical turning point for s- waves (now the 
angular momentum is not zero but rather h/2). 

We turn next to the eikonal approximation. This, may be considered as the high-energy 
limit of the WKB approximation. Expanding the WKB phase shift in V/E, and keeping the 
linear term, and further assuming r 2 = z 2 + b 2 , with b being the impact parameter, kb = A, 
we find for the eikonal Coulomb phase shift, a(b, k), the following, 

* m = -1/_>vpW < 62 > 

The integral in the above equation diverges at both extremes. Ref. [9] introduced a 
screening function that renders the integral finite. The screening function considered is 
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F(r) = @(a — r), where the screening length a is of very large, and taken to be a ^ b. Then, 



a(b,k) = --In 



which, with a ^> b, reduces to 

a(b, k) = 77 In 



b 

2a 



2a + 2Va 2 - b 2 
2a-2Va^¥ 



r/ In A — rjln(2ka) 



(63) 



(64) 



The above form is similar to Eq. (40), except for the change I — > A and the screening term, 
which depends only on energy and not on A. Glauber [U] considered also and exponential 
screening function of the form F(r) = exp —r/a, and found the following limiting expression 
for the Coulomb phase, 



a(b, k) = i] 



In 7 

2a ' 



(65) 



where 7 is Euler's constant, 7 = 0.57721.... Again, the main feature of Eq. (|40j) of 
a logarithmic dependence on b, and thus A is maintained. A Gaussian shaped screening 
function F(r) = exp [— r 2 /a 2 ] was considered by [TU], and the result found for the phase 
shift is, 

_ 1 
2a 2 



a(b, k) = r] 



In 



(66) 



Once again the main feature of Eq. (40) is maintained, namely, the logarithmic dependence 
on b. Clearly the discussion above demonstrates that the best route to follow to obtain 
well behaved and defined asymptotic forms for the Coulomb phase shift is to rely on the 
exact expression, Eq. 0, and use the uniform approximation method of Stirling's series to 
obtain the correct asymptotic form of the T function. 



IX. QUANTUM AND CLASSICAL COULOMB DEFLECTION FUNCTIONS 

An important theoretical entity which enters in any semi-classical treatment of scattering 
is the deflection function. In the continuous A = / + 1/2 limit, the deflection function is just 
the derivative of 2<r(A, k) with respect to A [I]. Thus, 



(67) 
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Using the WKB expression, Eq. (61), we obtain the well known classical Rutherford deflec- 
tion function, 

0(A) = 2tan~ 1 ^. (68) 
A 

The above should be compared to the exact "quantum" deflection function obtained from 
Ui of Eq. (|2]), which can be re- written as a recursion formula, viz, 

<Jl = ui-i + tan' 1 -. (69) 

Clearly one can define the derivative as merely the difference, a\— <Ji-\, and thus the quantum 
deflection function, Q q (l) is, 

e q (l) = 2[(7,-<t,_i] = 2tan- 1 ^. (70) 

The quantum and classical deflection functions agree if one makes the change I — > A, as 
expected. 



X. CONCLUSIONS 

In this paper we have revisited the literature on the point Coulomb phase shift, and 
sharpened the applicability and range of validity of several of the approximations usually 
employed to calculate it . In particular, the phase shift at zero angular momentum, 
<Jo(r]), as a function of the Sommerfeld parameter, 77, is calculated exactly and a simple 
analytic expression for it valid for small and large r] is found. This is useful when 
evaluating the phase shift as a function of angular momentum, I, using the exact formula, 
MV) = °o(v) + ELo tan" 1 (77/m), of Eq. §. 
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